Magnetic-field-dependent zero-bias diffusive anomaly in Pb oxide-n-InAs structures: 
Coexistence of two- and three-dimensional states 
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The results of experimental and theoretical studies of zero-bias anomaly (ZBA) in the Pb-oxide-n- 
InAs tunnel structures in magnetic field up to 6T are presented. A specific feature of the structures 
is a coexistence of the 2D and 3D states at the Fermi energy near the semiconductor surface. The 
dependence of the measured ZBA amplitude on the strength and orientation of the applied magnetic 
field is in agreement with the proposed theoretical model. According to this model, electrons tunnel 
into 2D states, and move diffusively in the 2D layer, whereas the main contribution to the screening 
comes from 3D electrons. 
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I. INTRODUCTION 

It is well known that the electron-electron interac- 
tion strongly influences the transport properties of dis- 
ordered conductors .El Even in the presence of weak dis- 
order (£ft~p 3> 1, where Ef is the Fermi energy, t p is 
the momentum relaxation time, and ft = 1) the electron- 
electron interaction suppresses the one-particle density of 
states at the Fermi level (diffusive anomaly). This leads 
to small deviations from Ohm's law in the current-voltage 
characteristics of a tunnel junction at small voltages V. 
The diffusive anomaly, which appears as a dip in the 
differential tunneling conductance G = dl/dV at zero 
bias, reveals itself in almost all tunneling experiments 
and has been studied in various tunneling structures .□ 
This should be distinguished from other nonlinearities 
of the current-voltage characteristics at low bias, which 
are due to different physical phenomena. The form of 
the diffusive zero-bias anomaly (ZBA) depends on the 
dimensionality: 5G(V) oc ln|V| for tunneling into two- 
dimensional (2D) conductors, and SG(V) oc y/\V\ for 
three-dimensional (3D) conductors. The width of the 
dip in the tunneling conductance is of order t^ 1 , and 
therefore cannot be observed in pure conductors. 

The first theoretical explanation of the diffusive ZBA, 
by Altshuler, Aronov and Lee in Refs. || and [|, was 
based on the diagrammatic perturbative method. For 
low-dimensional systems, this theory was subsequently 
extended beyond the perturbative treatment by Nazarov 
in Refs. || and |^ (see also Refs. ^,|[ where the realistic 
system is described by the coupling of the tunnel junc- 
tion with the effective electromagnetic environment), and 
later, by Levitov and Shytov in Rcf. ||. Nazarov also 
gave a transparent physical interpretation of the diffusive 



ZBA; immediately after an electron tunnels into the dif- 
fusive conductor and forms the distribution p(r, t), the 
system acquires an extra energy due to the interaction 
between this electron and the electrons in the conduc- 
tor (Coulomb barrier). Therefore, the electron density 
perturbation p(r, t) must spread under the Coulomb bar- 
rier in order to reach the final state. This process con- 
tributes a many-electron action S(t) (t ~ 1/eV is the 
time of spreading of the electron density perturbation) 
to the total tunneling action, and thereby, suppresses 
the tunneling current. In the regime of the Coulomb 
blockade effect (S(t) 3> 1) the tunneling current is al- 
most completely suppressed. Conversely, in good metals 
(sft p ^> 1) the Coulomb interaction is screened, so that 
the many-electron action is small, S(t) -C 1, and gives 
only small correction to the differential conductance. At 
T = 0, this takes the form 



1 dG 2e r /cm 1 

GdV = ~ Im \ S y UJ )U^-eV+iOi 



(1.1) 



The density of the tunneling electron /9 w (r) is given by 
a diffusion propagator (diffuson), whereas the electrody- 
namical potential W (r) which it excites is given by 



^(r) = / dr%(r,r> w (r'), 



(1.2) 



where V u (r, r') is the dynamically screened Coulomb, po- 
tential. The action S(ui) is then explicitly given byEl 



S(u) = - I (ir^(r)i(r). 



(1.3) 



This simple formula for the action displays an impor- 
tant role for the interface of the tunnel junction in the 
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ZBA in the case of tunneling into a 3D conductor. In- 
deed, after the electron tunnels through the barrier, it 
first appears on the surface of the conductor before prop- 
agating into the bulk. The surface of the conductor obvi- 
ously affects the spreading process of the electron density 
Pa,. Consequently, it affects the amplitude of the ZBA. 
For example, it can partially block the spreading of the 
electron into final state, giving rise to additional factor of 
2 in the amplitude of the anomalyEl (the electron prop- 
agates into the half space). This interface effect is even 
more pronounced in the presence of a magnetic field. The 
role of the magnetic field is twofold. It causes the Lorentz 
force, which blocks the spreading of the electron density, 
but it also induces a Hall voltage, which causes a drift 
along the interface, and thereby enhances the spread- 
ing. If the magnetic field B is perpendicular to the junc- 
tion interface, only the first effect contributes to the ZBA 
and gives a B 2 dependence of the ZBA.eI If the magnetic 
field is parallel to the junction interface, the two effects 
exactly cancel. This results in the strongly anisotropic 
magnetic field dependence of the ZBA predicted in Ref. 
pd| . Namely, the ZBA depends only on the component of 
the magnetic field perpendicular to the interface of the 
junction, as it would be in the case of tunneling into a 
2D conductor. This effect has probably been observed in 
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Motivated by this physical situation, we theoretically 
and experimentally investigated the ZBA in Pb-oxide- 
n-InAs structures in the presence of a magnetic field. 
We expected that the specific feature of these structures, 
namely, coexistence of 3D and 2D electron states near 
the surface of InAs, will strongly influence the ZBA and 
especially its magnetic field dependence. In particular, 
as the current in these structures can occur through the 
tunneling of electrons into both 2D and 3D states, the 
principle question which arises is whether the ZBA has 
2D or 3D character. The results of our study can be sum- 
marized as follows. The electrons tunnel into 2D states 
and move diffusively in a 2D layer, whereas the main 
contribution to the screening comes from 3D electrons. 
This gives rise to the unusual magnetic field dependence 
of the ZBA. When the magnetic field B is perpendicu- 
lar to the interface of the tunnel junction, the amplitude 
of the ZBA grows as B 2 in agreement with Ref. [|. The 
ZBA amplitude strongly depends on the orientation of 
the magnetic field, in agreement with Ref. |ll|. However, 
when the magnetic field lies in the plane of the junction 
interface, the magnetic field dependence does not disap- 
pear. Instead, the ZBA amplitude is linear in B. 



II. EXPERIMENTAL RESULTS 

The differential conductance G = dl/dV and its 
derivative (dG/dV) as a function of bias and magnetic 
field in Pb-oxide-n-InAs tunnel structures were investi- 
gated in a magnetic field up to 6 T at temperatures 4.2 
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FIG. 1. Bias dependencies of the differential conductance 
G for structure 1 at T= 1.6 K. 
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FIG. 2. Magnetic field dependence of dG/dV for different 
biases. The topmost curve is for Bin, and the others are 
for B || n (structure 1). All curves are at T= 4.2 K. 
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FIG. 3. Bias dependencies of (G(B)-G(0.1T))/G(V,0.1T) 
for different magnetic fields B || n at T= 4.2 K. 
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and 1.6 K. The tunnel structures were fabricated on n- 
InAs wafers with two different pairs of electron concen- 
tration and mobility: 9.7 x 10 17 cm -3 , and 1.5 x 10 4 



cm V s (structure 1); and 1.9 x 10 



17 



,-3 



and 



1.8 x 10 4 cm 2 V~ 1 s~ 1 (structure 2). Ultraviolet illumina- 
tion for 10-15 min in dry air was used to form the thin 
oxide, which served as a tunneling barrier. The Pb elec- 
trode was then evaporated through a mask. The tunnel 
contacts fabricated on each wafer were similar and re- 
sults are shown for one of several contacts fabricated on 
each wafer. The traditional modulation procedure was 
used for measuring the differential conductance and its 
derivative. Measurements showed that decreasing of the 
modulation amplitudes below 0.2 mV do not change the 
features in the G vs. V curves. Therefore, in all investi- 
gations the modulation amplitude was 0.2 mV. 

The dominant contribution to the current in the inves- 
tigated structure is a tunneling current. This is evident 
from the bias dependencies of the differential conduc- 
tance, which are shown on Fig. The structure of the 
curve for B = is the "superconducting anomaly" , con- 
nected with the superconducting gap in the one particle 
density of states in the metal electrode. At B > 0.06 T 
the superconductivity of Pb is destroyed and this struc- 
ture disappears completely. 

Oscillations in G and dG/dV as a function of V and B 
were observed for both B || n and Bin, where n is the 
normal to the plane of the tunnel junction (Figs. |^ and 
|3|). The tunneling conductance oscillations in such typ£S. 
of structures wete,|Gemprehensively studied in InAs,cj~t£l 
and in HgCdTe.Blia It was shown that in the structures 
based on InAs, an accumulation layer with 2D subbands 
exists near the barrier (Fig. [|). The tunneling conduc- 
tance is determined by tunneling into both 3D and 2D 
states of the semiconductor electrode. Having B || n 
leads to quantization of the spectrum of both 2D and 3D 
states. For this orientation of the magnetic field, the os- 
cillations in G are mainly due to the modulation of the 
density of 2D states. Having Bin does not quantize the 
energy spectrum of 2D states and the oscillations in G 
are only due to tunneling into 3D states. At fixed bias V, 
these oscillations are periodic in l/B. Therefore, using 
the Fourier transformation one can determine the funda- 
mental fields Bf and, consequently, the quasi- momenta 
k = ^/2eBf/cti of 2D and bulk states at the energy 
sf + eV. In addition, such data processing allows us to 
determine the energies of the bottoms of the conduction 
band and 2D subbands counted from the Fermi ener gy o f 
the semiconductor (for more details, see Rcfs.[ljand|l7|). 
Thus, we found that in structure 1 there are bulk states 
with ef — £ c = 115 meV and k 2 (£p) = 9.3 x 10 12 cm -2 , 
states of the ground 2D subband with ef — e° — 160 
meV and k 2 (sF) = 20.6 x 10 12 cmT 2 \ and states of 
the excited 2D subband with ef — e 1 — 120 meV and 



k 2 (£ F ) = 10.3 x 10 12 



rameters are Ef 



ef 



95 



cm . For the structure 2 these pa- 
- £c = 50 meV and k 2 (£ F ) = 3.1 x 10 12 , 
meV and fcg(e_p) = 7.6 x 10 12 , and 



Ef — e 1 ~ 55 meV and k 2 {£F) = 3.5 x 10 12 cm 2 . 

Now let us consider G vs. V curves in the vicin- 
ity of zero bias. The relative difference (G(V, B) — 
G(V, 0.1 T))/ G(V, 0.1 T) as a function of voltage for var- 
ious magnetic fields B || n is presented in Fig. ||. It is 
seen that increasing B gives rise to a dip in the con- 
ductance in the vicinity of V = 0, which is better seen 
when it falls between adjacent 2D Landau levels. This 
peculiarity is more pronounced in A (dG(V, B)/dV) = 
dG(V,B)/dV - dG(V,0.1 T)/dV vs. V curves (Fig. §a). 
To separate out the ZBA from the conductance oscilla- 
tions (due to Landau quantization), the following proce- 
dure was used; after taking the Fourier transformation 
(Fig. ||b) we cut out the components associated with the 
oscillations and then take the inverse Fourier transforma- 
tion (Fig. ||c). Such a procedure greatly helps in extract- 
ing the anomaly from the oscillations, but does not com- 
pletely separate the ZBA from the oscillations. There- 
fore, we cut out the part of the curve in the range ±5 mV 
in vicinity of V — (Fig. ||c), interpolate the rest of the 
curve by a smooth line, and then subtract this line from 
the initial curve shown in Fig. |^c. After integration, we 
obtain the ZBA in the tunneling conductance (Fig. ||d). 
(The correctness of such processing was verified by sepa- 
rating out the Gaussian shape from the simulating curve 
Ax sin(wi V + tpi) + A 2 s\n{uj 2 V + ip 2 ) + A 3 exp(-(U/A) 2 .) 

The magnetic field dependencies of the normalized am- 
plitude A = —SG/G\v=o and halfwidth of the ZBA are 
plotted in Fig. ^. It is seen that the halfwidth does not 
vary with the magnetic field within the experimental er- 
ror, whereas the amplitude of the ZBA significantly in- 
creases. The inset in Fig. ||a shows that the A vs. B 
dependence is close to i a B 2 . Similar results are ob- 
tained for the structure 2 (Fig. ^). In addition, one can 
see an oscillatory dependence of A on B, which appears 
at high magnetic fields. The minima of the oscillations 
are observed at those magnetic fields where the 2D Lan- 
dau levels cross the Fermi level. Thus, the origin of the 
oscillations of the ZBA amplitude is the Landau quan- 
tization of electron states in the 2D layer. The detailed 
investigation of this effect will be the subject of future 
work. Therefore, we will not concentrate on these oscil- 
lations in this paper. 

The angular dependence of the ZBA amplitude is plot- 
ted in Fig. U [if is the angle between B and n). One 
can see that the ZBA amplitude is strongly anisotropic. 
It drastically decreases when the magnetic field deviates 
from B || n, but it does not disappear at B _L n. The 
magnetic field dependence of the ZBA amplitude for this 
orientation is significantly weaker and is close to linear 
(Fig-0). 



3 



P b oxide 



n - I n A s 





i 

♦ e V 




1 1 1 1 1 1 II 

////// / / 


///// 






















FIG. 4. Energy diagram of the Pb-oxide-n-InAs tunnel 
structure. e° and e 1 are the energies of the bottom of the 
ground and excited 2D subbands respectively, and e c is the 
energy of the bottom of the conduction band. 
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FIG. 5. (a) The bias dependence of A(dG(V, B)/dV) = 
dG{V,B)/dV - dG(V, 0.1 T)/dV for B = 5.3 T and B || n. 
(b) Fourier transform of the upper curve. The bar shows the 
region which was cut out. (c) Result of inverse Fourier trans- 
form. The dotted curve is the interpolation after the central 
region was cut out. (d) Reconstructed zero-bias anomaly after 
the processing described in the text. 
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FIG. 6. Magnetic field dependencies of the amplitude (a) 
and halfwidth (b) of the ZBA, for B || n. The inset shows the 
A vs. B 2 dependence. All points are at T= 4.2 K. 




FIG. 7. Magnetic field dependencies of the ZBA amplitude 
for B || n (open squares) and B _!_ n (full squares) for struc- 
ture 2 at T= 4.2 K. The arrows indicate the magnetic fields, 
for which the 2D Landau levels coincide with the Fermi level. 



III. DISCUSSION 

The energy scale of the diffusive anomaly caused by 
electron-electron interaction in dirty conductors is t^ 1 . 
In our structures, t" 1 is estimated from the mobility to 
be about 2 meV, whereas the halfwidth of the ZBA is 1 
meV (Fig. |b). Thus, we suppose that the ZBA observed 
in our experiment is just the diffusive anomaly. The spe- 
cific feature of the investigated structures is the coexis- 
tence of 2D and 3D electrons near the barrier. Therefore 
the basic question is whether the zero-bias anomaly is 
due to the interaction of 2D or 3D electrons. In princi- 
ple, the ZBA has different form for tunneling into 2D a nd 
3D states: SG(V) oc ln\V\ for 2D and 5G(V) oc y/\V~\ for 
3D. However, the comparison of fits to experimental data 
in Fig. [)]does not allow us to distinguish between the two 
forms of the ZBA. 

On one hand, the main part of the tunneling conduc- 
tance is due to tunneling into the 2D states. This follows 
from the theoretical calculation of the tunneling conduc- 
tance for the investigated structures carried ouj_in the 
framework of the transfer Hamiltonian method.Efl Such 
a calculation shows that the tunneling conductance due 
to tunneling into the 2D states is larger by about a fac- 
tor of 5 than that for tunneling into 3D states. This 
conclusion is also supported by the fact that the ampli- 
tude of oscillations of the tunneling conductance caused 
by the Landau quantization is significantly larger in the 
case of tunneling into 2D states (at B || n) than in the 
case of tunneling into 3D states (B _L n) (Fig. ||). In 
addition, the ZBA amplitude has typical for 2D systems 
strong dependence on the magnetic field orientation, i.e. 
it is determined mainly by the normal component of the 
magnetic field (Fig. |J). Therefore, one can surmise that 
the ZBA has 2D character. 

On the other hand, the strong angular dependence of 
the ZBA is ambiguous evidence of the 2D nature of the 
ZBA. Indeed, in Ref. |ll| it was demonstrated that for 
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B ll n <P (degrees) Bin 

FIG. 8. Angular dependence of the ZBA amplitude, for 
B=5.5 T. The crosses and squares are data for structures 1 
and 2 respectively. The dashed curve is the result of the 
calculation described in the text. The solid curve corresponds 
to the 3D case. 

tunneling into 3D states the amplitude of the ZBA is 
given by A(B) oc 1 +w^Tp cos 2 (cp) (where lo c is cyclotron 
frequency), i.e. it depends only on the component of the 
magnetic field perpendicular to the interface of the junc- 
tion, as it would be in the case of tunneling into a 2D 
conductor. In the limit w^r 2 ^> 1 this leads to the 
strong angular dependence of the ZBA amplitude. Al- 
though the strong anisotropy of the ZBA is observed in 
our experiment for uj^t 2 ~ 50 at B = 5.5 T, the curve 
1 + u>c T p cos 2 (<p) does not fit well with the experimental 
data (see the Fig. ||). Moreover, the ZBA amplitude is 
linear in B when the magnetic field lies in the plane of 
the junction interface (see Fig. |J). 

Thus, the magnetic field dependence of the ZBA in the 
investigated structures does not completely agree with 
either the 3D or 2D nature of the ZBA. We would like 
to stress however, that our experimental set up is not 
usual for studying the diffusive ZBA. Traditionally, the 
2D metallic layer in tunnel junctions is electrically iso- 
lated from the 3D electrode (or another 2D layer). In this 
case the charge relaxation is two-dimensional and the in- 
teraction is partially screened by the 3D metal, so that 
its strength is defined by the distance A between 2D and 
3D electrodes. The correction to the differential conduc- 
tance then has the forrriiJ SG(V) ~ ln(aA/r|j) ln(eW p ), 
where a is the width of the 2D layer, and ru is the Debye 
radius. In addition, it is assumed that this formula holds 
for aA/r|, 3> 1 and A/t\d ^> 1. Thus, there are two rea- 
sons which make our experimental set up different from 
the usual one, and the above formula nonapplicable to 
our case. The specific feature of the investigated struc- 
tures is coexistence of 3D and 2D electron states near the 
surface of the semiconductor (see Fig. 0). Thus, formally 
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FIG. 9. Tunneling conductance near V = versus In |V| 
(a) and versus | V"| 5 (b) at T= 1.6 K, B=5.5 T The two sets 
of data points correspond to different signs of the bias. 

in our case A = 0. Second, and this is most important, 
in our experiment the 2D electron system and 3D metal 
are not electrically isolated. 

In the next section we show that the following sce- 
nario of tunneling is realized in this structures. After 
tunneling, the electron moves diffusively in the 2D layer 
and forms the 2D distribution p w (r) at the surface of the 
semiconductor. It immediately pushes other electrons 
into the bulk, so that the total charge becomes zero for a 
short time of order of the inverse plasma frequency. Thus, 
the relaxation of the total charge takes three-dimensional 
form, and this should lead to y/V -dependence of the dif- 
ferential conductance usual for the 3D ZBA. The dimen- 
sionality of p u (r) does not affect the voltage dependence 
of the differential conductance. However, it leads to the 
unusual magnetic field dependence of the ZBA discussed 
above. 



IV. THEORETICAL MODEL AND 
COMPARISON WITH EXPERIMENTAL DATA 

In our theoretical analysis we make two assumptions. 
First, we assume that the probability of electrons tunnel- 
ing into 2D states near the surface of the semiconductor 
is much greater than the probability of tunneling into 
3D bulk states. This follows from the analysis of the ex- 
perimental results in the previous sections. The second 
assumption is that the tunneling electron (though being 
screened by the other electrons of the system) remains 
in the 2D well for a physically relevant time, i.e. for the 
time t ~ 1/eV (see discussion in Sec. [|) before escaping 
into the bulk. This last time ranges from t v to (fcT) -1 
in the ZBA regime, and thus in our experiment it may 
exceed the momentum relaxation time by factor of 10. 
On the other hand, the escape of electrons from the 2D 
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well to the bulk is determined by the ionized impurity 
scattering, which is main scattering mechanism at low 
temperatures. This mechanism is strongly anisotropic - 
the small angle scattering dominates. Together with re- 
quirement of large momentum transfer for 2D— »3D tran- 
sition this leads to the fact that the 2D — > 2D transition 
rate W22 ~ ( t p) _1 is larger than the 2D — ► 3D tran- 
sition rate, W23. The calculations with wave functions 
and screening radius corresponding to the investigated 
structures carried out in the same manner as in Ref. |20| 
gives W22/W23 ~ 15. Thus, the second assumption is 
justified. 

The width of the 2D well is of order Ap and is much 
smaller than the mean free path I3 in the bulk of the 
conductor. Therefore, the 2D well can be thought of 
as a 5-layer with respect to the physically relevant length 
scale. This means that after tunneling the electron forms 
a 2D density distribution p w (r) = Q UJ (R)S(z) localized on 
the surface z = of the conductor (here, R = (x, y) is the 
coordinate on the surface, and r = (R, z)). Then, due 
to the interaction of this electron with the ones forming 
both 2D and 3D liquids, the electrodynamic potential <f> u 
is excited. 



Instead of a direct calculation of the integral (1.2) for 
(f> u , we follow Ref. [ll] and use the electroneutrality prin- 
ciple. We assume that the density of the tunneling elec- 
tron is completely screened on the distance of the or- 
der of Debye radius t\d, so that the induced charge den- 
sity is p u (t) = -p w (rL== -Q UJ (R)S(z). Taking into ac- 
count charging effectsEj gives only corrections of order 
r 2 D /1VI3 -C 1 (w is the thickness of the tunneling barrier) , 
which we neglect here. We also assume that the Pb elec- 
trode, being a good metal, does not contribute to the 
action S(uj). Therefore, after Fourier transformation the 



integral (1.3) can be represented in the following form: 



1 

8^2 



dkQ_ w (-k)$ w (k), 



(4.1) 



where $ w (k) = w (k, z)\ z=0 . 

The density Q u obeys the 2D diffusion equation in 
imaginary time 

\u\Q u - £> 2 V^Q W = -e sign(io)5(K - R ), (4.2) 

with the diffusion coefficient D 2 = D 2 (B), which de- 
pends only on the z-component of the magnetic field 
B = (B,<p): 



D 2 (B) 



D (0) 
^2 



1 + CJ 2 t| cos 2 <p 



(4.3) 



Here, w c is the cyclotron frequency and D { 2 0) = D 2 (0). We 
have introduced the new notation, t 2 , for the momentum 
relaxation time of 2D electrons to distinguish it from that 
of 3D electrons. After Fourier transformation, equation 
(4.2) can be immediately solved, 



Q«(k) = - 



e signitu) 



(4.4) 



To calculate the potential <& w we formulate and then 
solve the system of equations for the dynamics of the 
induced charge density p u — —Q ul 5{z). This dynam- 
ics is controlled by the transport along the 2D layer, as 
well as by the nonzero current perpendicular to the layer 
j n (oj, R) . The conservation of charge (in imaginary time) 
reads, 



MQa 



D 2 V 2 R Q U 



J 1 



(4.5) 



where we introduced the 2D conductivity (72(B) = 
e 2 ^2-D2(B), and v 2 is the Fermi density of 2D states. 
On the left hand side of this equation the second and 
third terms are the divergences of the diffusion and elec- 
trical currents respectively. The first two terms of this 
equation coincide with the left-hand side of Eq. (|4.2| ) for 
the diffusion propagator Q u . This is precisely the reason 
for the cancellation of the diffusion pole discussed below. 
Now, we can use this fact to eliminate Q u from the last 
equation: 



Jr. 



o^V 2 ^^ = — e sign(u))5(H — Rq)- 



(4.6) 



On the other hand, in the bulk of the conductor the 
charge is not accumulated, Pu,\ z >o — 0. The conservation 
of charge then leads to Vj(w,r) = and n-j(u>, r)| z= o = 
j n {^, R), where j is the density of current in the bulk of 
the conductor. These two equations can be expressed in 
terms of d>^ : 



V r -<rV r w (R, z) = 0, 



n-<7V r <^(R, z)\ z =o = -in(w,R), 



(4.7) 



(4.8) 



where a = (7(B) is the conductivity tensor in the bulk 
of the conductor. If the magnetic field is perpendicular 
to the surface of the conductor (ip = 0), the conductivity 
tensor takes the simple form: 



'vv 



2J2 ■ 



■W'T. 



(0) 



c'3 



Or 



-sign(u>] 



(0) 



(4.9) 
(4.10) 



and all other elements vanish. Here, is the con- 
ductivity of 3D electrons in the case of zero magnetic 
field, and t 3 is the momentum relaxation time of 3D elec- 
trons. In the case of arbitrary magnetic field orientation, 
a can be calculated by the rotation over the angle ip, 
a(B) = U(-v)a(B)L= U(<p). 

Solving Eqs. (4.6,4.7,4.8) simultaneously, we get, 



Mk) = - 



e sign(u>) 
(73 0) Z + a2k 2 ' 



where the function Z(B,k) is given by 



Z(B,k) 



A: 2 cos 2 a + ky + ik x sin a 



(4.11) 



(4.12) 
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with the angle a defined by 



sin a 



smtp 



sign(u>) 



(4.13) 



Finally, we s ubstitute Q u from (4.4) and <& w from (4.11) 
into Eq. (4.1), and arrive at the following formula for the 
action: 



' SV) ' 8^ 2 ./ (M Lhk-)(n-;;; z-:- ^k-) 



dk 



(4.14) 



We would like to mention th e can cellati on o f the diffu- 
sion pol e in the expressions (4.11) and (4.14) (compare 
to Eqs. ( |l.2|jl.3| )). This well known fact (see for example 
Ref. [l]) physically means that after tunneling, the accom- 
modation of the charge is entirely governed by plasmon 
modes. 



In the equation (4.14) two terms in the denominator 



of the integrand, <7 3 Z and 02 k 2 , are contributions from 
the screening of the tunneling electron by 3D and 2D 
electrons, respectively. In principle, one can expect to 
observe the crossover from 3D ZBA with S(uj) ~ 1/V^ 
(g = <7 ( ° } Z/a 2 k 2 > 1) to 2D ZBA with S(u) ~ 1/w 
(<7 <C 1). In our case, however, the s creen ing by 2D 
electrons is weak. Indeed, the integral fl4.14| ) should be 
evaluated for ui ~ I/V2. For p = 0, the simple estimate 
then gives g ~ kph 3> 1 in the case of zero magnetic field, 
and g ~ k F l 2 ^ 1 in strong magnetic field, lo c t 2 ^ ^> 1 
(here, l 2 is the mean free path in the 2D layer). For 
if = n/2 we have g ~ Ef /u c ^> lin strong magnetic field. 
Therefore, we can neglect screening by 2D electrons. This 
equally means that in Eq. ( j4.E| ), for the induced charge 
density Q Ul we can neglect the 2D current, -ff 2 Vn$ u , 
compared to the 3D current, j n (given by Eq. (4.8)). 



Thus, after tunneling the charge relaxation process has 
3D character. This leads to a y/V dependence of the 
differential conductance, as is usual for 3D. We show this 
next. 

Omitting t he ter m (72 k 2 in the denominator in the right 
hand side of (4.14), we carry out the integration over k 
and obtain, 



The correction to the differential conductance then takes 
the form SG ~ W. 

Next we concentrate on the magnetic field dependence 
of the ZBA. The fact that the magnetic field and V de- 
pendencies of the differential conductance are completely 
factorized allows us to represent the normalized ampli- 
tude of the ZBA, A(B,ip) = —5G /G\y =o, in a simple 
form. To do this, we integrate Eq. ( 4.17 ) over V and cut 
the integral at eV ~ ', 
result: 



We then arrive at the following 



A(B,<p)=A F(B,<p) i (4.18) 

where Ao = nXp/^h, and k is a dimcnsionless num- 
ber of order 1. The equation (4.18), together with Eq. 



(4.16), represents the general result which is valid for an 
arbitrary magnetic field. Now we consider the most in- 
teresting case of strong magnetic field, w c T3 3> 1. In this 
limit sin 2 a = sin 2 ip and introducing the dimensionless 
parameter h — lo c t 2 we can write, 



A(B,(p) = k (Xp/hf h\fl + h 2 cos 2 ^(sin ip). (4.19) 

We are now in a position to compare the result of our 
theoretical analysis with the experimental data. When 
the magnetic field is perpendicular to the interface of the 
tunnel junction, ip — 0, we have (co c T2 S> 1): 



T ( R~ c 



(4.20) 



where R c is the cyclotron radius. The amplitude of the 
anomaly thus goes as B 2 . As it is clearly seen from Figs. 
^ and (7], the experimental data for B |j n are close to this 
dependence for both structures. Remarkably, the ampli- 
tude ofthe anomaly does not depend on T2 and T3 (see 
Eq. ( 4.20) )). However, contrary to one's first expectation, 
the ZBA cannot be observed in a perfect 3D metal. Al- 
though the amplitude of the anomaly stays constant with 
7"2,3 ~~ * 00) the dip of the tunneling conductance gets nar- 
rower (its width is given by t 2 x ) and finally shrinks. 

When the magnetic field is parallel to the junction in- 
terface, ip — 7r/2, from Eq. (4.1£) we obtain, 



--F(B,<p), 



(4.15) 



F{B, <p) = sj (1 + lu 2 c t 2 ) (1 + u 2 t 2 cos 2 ^)E(sm a), 



(4.16) 



tt/2 



where E(s) = j dOy/l - s 2 sin 2 is the complete ellip- 



tic integr al. S ubstituting the action S(u>) from Eq. (4.15) 
into Eq. (1.1), we arrive at the final result, 



1 dG 
GdV 



2ir 2 a ( °'>j2eVD 2 0) 



-F(B, tp). 



(4.17) 



A(B,ir/2) 



h 



(4.21) 



i.e. the amplitude of ZBA is a linear function of magnetic 
field. The same dependence is observed experimentally 
(see Fig. 0). 

Finally, we can keep the amplitude of the magnetic 
field constant and study the angular dependence of the 
ZBA. From Eq. ([4.190 it follows that 



A(B,0) 



7m 



h 2 cos 2 ipE(smip). 



(4.22) 



This dependence rSjid the dependence corresponding to 
the pure 3D caseEJ are p lotted in Fig. ||. One can see 
that the expression ( |4. 22| ) is in excellent agreement with 
the experimental data without any fitting parameters. 



7 



V. CONCLUSION 

We have presented the results of experimental and the- 
oretical studies of the zero-bias anomaly (ZBA) in tunnel 
structures with 2D and 3D electron states coexisting near 
the semiconductor surface. It has been shown that the 
specific scenario of tunneling realized in this structures 
is: (i) electrons tunnel mainly into 2D states, (ii) imme- 
diately after tunneling, the electrons move diffusively in 
a 2D layer, and (iii) the main contribution to screening 
comes from the 3D electrons and, as a result, the charge 
relaxation has a 3D character. 

This leads to the peculiar features of the magnetic 
field dependence of the ZBA amplitude. When the mag- 
netic field is perpendicular to the interface of the tunnel 
junction, the ZBA amplitude grows as B 2 , in agreement 
with Ref. |i[ Although the magnetic field dependence has 
strong anisotropy, as predicted in Ref. |ll|, it does not 
disappear completely when the magnetic field lies in the 
plane of the junction interface. Instead, the ZBA ampli- 
tude is linear in B. 

The experimental data show that the ZBA amplitude 
oscillates with the magnetic field. The origin of the os- 
cillations is the Landau quantization of electron states in 
the 2D layer. The detailed investigation of this effect will 
be the subject of future work. 
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